Genetic incompatibilities in reciprocal hybrids between populations of Tigriopus californicus with low to moderate mitochondrial sequence divergence.

All mitochondrial-encoded proteins and RNAs function through interactions with nuclear-encoded proteins, which are critical for mitochondrial performance and eukaryotic fitness. Coevolution maintains inter-genomic (i.e., mitonuclear) compatibility within a taxon, but hybridization can disrupt coevolved interactions, resulting in hybrid breakdown. Thus, mitonuclear incompatibilities may be important mechanisms underlying reproductive isolation and, potentially, speciation. Here we utilize Pool-seq to assess the effects of mitochondrial genotype on nuclear allele frequencies in fast- and slow-developing reciprocal inter-population F2 hybrids between relatively low-divergence populations of the intertidal copepod Tigriopus californicus. We show that mitonuclear interactions lead to elevated frequencies of coevolved (i.e., maternal) nuclear alleles on two chromosomes in crosses between populations with 1.5% or 9.6% fixed differences in mitochondrial DNA nucleotide sequence. However, we also find evidence of excess mismatched (i.e., noncoevolved) alleles on three or four chromosomes per cross, respectively, and of allele frequency differences consistent with effects involving only nuclear loci (i.e., unaffected by mitochondrial genotype). Thus, our results for low-divergence crosses suggest an underlying role for mitonuclear interactions in variation in hybrid developmental rate, but despite substantial effects of mitonuclear coevolution on individual chromosomes, no clear bias favoring coevolved interactions overall.

Developmental rate displays effects of inheritance but not of sex in interpopulation hybrids of Tigriopus californicus.

Coevolved genetic interactions within populations can be disrupted by hybridization resulting in loss of fitness in hybrid individuals (i.e., hybrid breakdown). However, the extent to which variation in fitness-related traits among hybrids is inherited across generations remains unclear, and variation in these traits may be sex-specific in hybrids due to differential effects of genetic incompatibilities in females and males. Here we present two experiments investigating variation in developmental rate among reciprocal interpopulation hybrids of the intertidal copepod Tigriopus californicus. Developmental rate is a fitness-related trait in this species that is affected by interactions between mitochondrial-encoded and nuclear-encoded genes in hybrids that result in variation in mitochondrial ATP synthesis capacities. First, we show that F2 -hybrid developmental rate is equivalent in two reciprocal crosses and is unaffected by sex, suggesting that breakdown of developmental rate is likely experienced equally by females and males. Second, we demonstrate that variation in developmental rate among F3 hybrids is heritable; times to copepodid metamorphosis of F4 offspring of fast-developing F3 parents (12.25 ± 0.05 days, µ ± SEM) were significantly faster than those of F4 offspring of slow-developing parents (14.58 ± 0.05 days). Third, we find that ATP synthesis rates in these F4 hybrids are unaffected by the developmental rates of their parents, but that mitochondria from females synthesize ATP at faster rates than mitochondria from males. Taken together, these results suggest that sex-specific effects vary among fitness-related traits in these hybrids, and that effects likely associated with hybrid breakdown display substantial inheritance across hybrid generations.

Differential gene expression and mitonuclear incompatibilities in fast- and slow-developing interpopulation Tigriopus californicus hybrids.

Mitochondrial functions are intimately reliant on proteins and RNAs encoded in both the nuclear and mitochondrial genomes, leading to inter-genomic coevolution within taxa. Hybridization can break apart coevolved mitonuclear genotypes, resulting in decreased mitochondrial performance and reduced fitness. This hybrid breakdown is an important component of outbreeding depression and early-stage reproductive isolation. However, the mechanisms contributing to mitonuclear interactions remain poorly resolved. Here, we scored variation in developmental rate (a proxy for fitness) among reciprocal F2 interpopulation hybrids of the intertidal copepod Tigriopus californicus and used RNA sequencing to assess differences in gene expression between fast- and slow-developing hybrids. In total, differences in expression associated with developmental rate were detected for 2925 genes, whereas only 135 genes were differentially expressed as a result of differences in mitochondrial genotype. Upregulated expression in fast developers was enriched for genes involved in chitin-based cuticle development, oxidation-reduction processes, hydrogen peroxide catabolic processes and mitochondrial respiratory chain complex I. In contrast, upregulation in slow developers was enriched for DNA replication, cell division, DNA damage and DNA repair. Eighty-four nuclear-encoded mitochondrial genes were differentially expressed between fast- and slow-developing copepods, including 12 subunits of the electron transport system (ETS) which all had higher expression in fast developers than in slow developers. Nine of these genes were subunits of ETS complex I. Our results emphasize the major roles that mitonuclear interactions within the ETS, particularly in complex I, play in hybrid breakdown, and resolve strong candidate genes for involvement in mitonuclear interactions.

Loss of mitochondrial performance at high temperatures is correlated with upper thermal tolerance among populations of an intertidal copepod.

Environmental temperatures have pervasive effects on the performance and tolerance of ectothermic organisms, and thermal tolerance limits likely play key roles underlying biogeographic ranges and responses to environmental change. Mitochondria are central to metabolic processes in eukaryotic cells, and these metabolic functions are thermally sensitive; however, potential relationships between mitochondrial function, thermal tolerance limits and local thermal adaptation in general remain unresolved. Loss of ATP synthesis capacity at high temperatures has recently been suggested as a mechanistic link between mitochondrial function and upper thermal tolerance limits. Here we use a common-garden experiment with seven locally adapted populations of intertidal copepods (Tigriopus californicus), spanning approximately 21.5° latitude, to assess genetically based variation in the thermal performance curves of maximal ATP synthesis rates in isolated mitochondria. These thermal performance curves displayed substantial variation among populations with higher ATP synthesis rates at lower temperatures (20-25 °C) in northern populations than in southern populations. In contrast, mitochondria from southern populations maintained ATP synthesis rates at higher temperatures than the temperatures that caused loss of ATP synthesis capacity in mitochondria from northern populations. Additionally, there was a tight correlation between the thermal limits of ATP synthesis and previously determined variation in upper thermal tolerance limits among populations. This suggests that mitochondria may play an important role in latitudinal thermal adaptation in T. californicus, and supports the hypothesis that loss of mitochondrial performance at high temperatures is linked to whole-organism thermal tolerance limits in this ectotherm.

Strong selective effects of mitochondrial DNA on the nuclear genome.

Oxidative phosphorylation, the primary source of cellular energy in eukaryotes, requires gene products encoded in both the nuclear and mitochondrial genomes. As a result, functional integration between the genomes is essential for efficient adenosine triphosphate (ATP) generation. Although within populations this integration is presumably maintained by coevolution, the importance of mitonuclear coevolution in key biological processes such as speciation and mitochondrial disease has been questioned. In this study, we crossed populations of the intertidal copepod Tigriopus californicus to disrupt putatively coevolved mitonuclear genotypes in reciprocal F2 hybrids. We utilized interindividual variation in developmental rate among these hybrids as a proxy for fitness to assess the strength of selection imposed on the nuclear genome by alternate mitochondrial genotypes. Developmental rate varied among hybrid individuals, and in vitro ATP synthesis rates of mitochondria isolated from high-fitness hybrids were approximately two-fold greater than those of mitochondria isolated from low-fitness individuals. We then used Pool-seq to compare nuclear allele frequencies for high- or low-fitness hybrids. Significant biases for maternal alleles were detected on 5 (of 12) chromosomes in high-fitness individuals of both reciprocal crosses, whereas maternal biases were largely absent in low-fitness individuals. Therefore, the most fit hybrids were those with nuclear alleles that matched their mitochondrial genotype on these chromosomes, suggesting that mitonuclear effects underlie individual-level variation in developmental rate and that intergenomic compatibility is critical for high fitness. We conclude that mitonuclear interactions can have profound impacts on both physiological performance and the evolutionary trajectory of the nuclear genome.

Variation in developmental temperature alters adulthood plasticity of thermal tolerance in Tigriopus californicus.

In response to environmental change, organisms rely on both genetic adaptation and phenotypic plasticity to adjust key traits that are necessary for survival and reproduction. Given the accelerating rate of climate change, plasticity may be particularly important. For organisms in warming aquatic habitats, upper thermal tolerance is likely to be a key trait, and many organisms express plasticity in this trait in response to developmental or adulthood temperatures. Although plasticity at one life stage may influence plasticity at another life stage, relatively little is known about this possibility for thermal tolerance. Here, we used locally adapted populations of the copepod Tigriopus californicus to investigate these potential effects in an intertidal ectotherm. We found that low latitude populations had greater critical thermal maxima (CTmax) than high latitude populations, and variation in developmental temperature altered CTmax plasticity in adults. After development at 25°C, CTmax was plastic in adults, whereas no adulthood plasticity in this trait was observed after development at 20°C. This pattern was identical across four populations, suggesting that local thermal adaptation has not shaped this effect among these populations. Differences in the capacities to maintain ATP synthesis rates and to induce heat shock proteins at high temperatures, two likely mechanisms of local adaptation in this species, were consistent with changes in CTmax owing to phenotypic plasticity, which suggests that there is likely mechanistic overlap between the effects of plasticity and adaptation. Together, these results indicate that developmental effects may have substantial impacts on upper thermal tolerance plasticity in adult ectotherms.

Mitochondria, sex and variation in routine metabolic rate.

Variation in the metabolic costs associated with organismal maintenance may play a key role in determining fitness, and thus these differences among individuals are likely to be subject to natural selection. Although the evolvability of maintenance metabolism depends on its underlying genetic architecture, relatively little is known about the nature of genetic variation that underlies this trait. To address this, we measured variation in routine metabolic rate (MO2routine ), an index of maintenance metabolism, within and among three populations of Atlantic killifish, Fundulus heteroclitus, including a population from a region of genetic admixture between two subspecies. Polygenic association tests among individuals from the admixed population identified 54 single nucleotide polymorphisms (SNPs) that were associated with MO2routine , and these SNPs accounted for 43% of interindividual variation in this trait. However, genetic associations with MO2routine involved different SNPs if females and males were analysed separately, and there was a sex-dependent effect of mitochondrial genotype on variation in routine metabolism. These results imply that there are sex-specific genetic mechanisms, and potential mitonuclear interactions, that underlie variation in MO2routine . Additionally, there was evidence for epistatic interactions between 17% of the possible pairs of trait-associated SNPs, suggesting that epistatic effects on MO2routine are common. These data demonstrate not only that phenotypic variation in this ecologically important trait has a polygenic basis with considerable epistasis among loci, but also that these underlying genetic mechanisms, and particularly the role of mitochondrial genotype, may be sex-specific.

Mitochondrial Ecophysiology: Assessing the Evolutionary Forces That Shape Mitochondrial Variation.

The mitonuclear species concept hypothesizes that incompatibilities between interacting gene products of the nuclear and mitochondrial genomes are a major factor establishing and maintaining species boundaries. However, most of the data available to test this concept come from studies of genetic variation in mitochondrial DNA, and clines in the mitochondrial genome across contact zones can be produced by a variety of forces. Here, we show that using a combination of population genomic analyses of the nuclear and mitochondrial genomes and studies of mitochondrial function can provide insight into the relative roles of neutral processes, adaptive evolution, and mitonuclear incompatibility in establishing and maintaining mitochondrial clines, using Atlantic killifish (Fundulus heteroclitus) as a case study. There is strong evidence for a role of secondary contact following the last glaciation in shaping a steep mitochondrial cline across a contact zone between northern and southern subspecies of killifish, but there is also evidence for a role of adaptive evolution in driving differentiation between the subspecies in a variety of traits from the level of the whole organism to the level of mitochondrial function. In addition, studies are beginning to address the potential for mitonuclear incompatibilities in admixed populations. However, population genomic studies have failed to detect evidence for a strong and pervasive influence of mitonuclear incompatibilities, and we suggest that polygenic selection may be responsible for the complex patterns observed. This case study demonstrates that multiple forces can act together in shaping mitochondrial clines, and illustrates the challenge of disentangling their relative roles.

Variation in Thermal Tolerance and Its Relationship to Mitochondrial Function Across Populations of Tigriopus californicus.

Variation in thermal tolerance plays a key role in determining the biogeographic distribution of organisms. Consequently, identifying the mechanistic basis for thermal tolerance is necessary for understanding not only current species range limits but also the capacity for range limits to shift in response to climate change. Although variation in mitochondrial function likely contributes to variation in thermal tolerance, the extent to which mitochondrial function underlies local thermal adaptation is not fully understood. In the current study, we examine variation in thermal tolerance and mitochondrial function among three populations of the intertidal copepod Tigriopus californicus found across a latitudinal thermal gradient along the coast of California, USA. We tested (1) acute thermal tolerance using survivorship and knockdown assays, (2) chronic thermal tolerance using survivorship of nauplii and developmental rate, and (3) mitochondrial performance at a range of temperatures using ATP synthesis fueled by complexes I, II, and I&II, as well as respiration of permeabilized fibers. We find evidence for latitudinal thermal adaptation: the southernmost San Diego population outperforms the northernmost Santa Cruz in measures of survivorship, knockdown temperature, and ATP synthesis rates during acute thermal exposures. However, under a chronic thermal regime, survivorship and developmental rate are more similar in the southernmost and northernmost population than in the mid-range population (Abalone Cove). Though this pattern is unexpected, it aligns well with population-specific rates of ATP synthesis at these chronic temperatures. Combined with the tight correlation of ATP synthesis decline and knockdown temperature, these data suggest a role for mitochondria in setting thermal range limits and indicate that divergence in mitochondrial function is likely a component of adaptation across latitudinal thermal gradients.

Patterns of alternative splicing in response to cold acclimation in fish.

Phenotypic plasticity is an important aspect of an organism's response to environmental change that often requires the modulation of gene expression. These changes in gene expression can be quantitative, as a result of increases or decreases in the amounts of specific transcripts, or qualitative, as a result of the expression of alternative transcripts from the same gene (e.g. via alternative splicing of pre-mRNAs). Although the role of quantitative changes in gene expression in phenotypic plasticity is well known, relatively few studies have examined the role of qualitative changes. Here, we use skeletal muscle RNA-seq data from Atlantic killifish (Fundulus heteroclitus), threespine stickleback (Gasterosteus aculeatus) and zebrafish (Danio rerio) to investigate the extent of qualitative changes in gene expression in response to cold acclimation. Fewer genes demonstrated alternative splicing than differential expression as a result of cold acclimation; however, differences in splicing were detected for 426 to 866 genes depending on species, indicating that large numbers of qualitative changes in gene expression are associated with cold acclimation. Many of these alternatively spliced genes were also differentially expressed, and there was functional enrichment for involvement in muscle contraction among the genes demonstrating qualitative changes in response to cold acclimation. Additionally, there was a common group of 29 genes with cold-acclimation-mediated changes in splicing in all three species, suggesting that there may be a set of genes with expression patterns that respond qualitatively to prolonged exposure to cold temperatures across fishes.

Integrative Population and Physiological Genomics Reveals Mechanisms of Adaptation in Killifish.

Adaptive divergence between marine and freshwater (FW) environments is important in generating phyletic diversity within fishes, but the genetic basis of this process remains poorly understood. Genome selection scans can identify adaptive loci, but incomplete knowledge of genotype-phenotype connections makes interpreting their significance difficult. In contrast, association mapping (genome-wide association mapping [GWAS], random forest [RF] analyses) links genotype to phenotype, but offer limited insight into the evolutionary forces shaping variation. Here, we combined GWAS, RF, and selection scans to identify loci important in adaptation to FW environments. We utilized FW-native and brackish water (BW)-native populations of Atlantic killifish (Fundulus heteroclitus) as well as a naturally admixed population between the two. We measured morphology and multiple physiological traits that differ between populations and may contribute to osmotic adaptation (salinity tolerance, hypoxia tolerance, metabolic rate, body shape) and used a reduced representation approach for genome-wide genotyping. Our results show patterns of population divergence in physiological capabilities that are consistent with local adaptation. Population genomic scans between BW-native and FW-native populations identified genomic regions evolving by natural selection, whereas association mapping revealed loci that contribute to variation for each trait. There was substantial overlap in the genomic regions putatively under selection and loci associated with phenotypic traits, particularly for salinity tolerance, suggesting that these regions and genes are important for adaptive divergence between BW and FW environments. Together, these data provide insight into the mechanisms that enable diversification of fishes across osmotic boundaries.

Tolerance traits related to climate change resilience are independent and polygenic.

The resilience of organisms to climate change through adaptive evolution is dependent on the extent of genetically based variation in key phenotypic traits and the nature of genetic associations between them. For aquatic animals, upper thermal tolerance and hypoxia tolerance are likely to be a important determinants of sensitivity to climate change. To determine the genetic basis of these traits and to detect associations between them, we compared naturally occurring populations of two subspecies of Atlantic killifish, Fundulus heteroclitus, that differ in both thermal and hypoxia tolerance. Multilocus association mapping demonstrated that 47 and 35 single nucleotide polymorphisms (SNPs) explained 43.4% and 51.9% of variation in thermal and hypoxia tolerance, respectively, suggesting that genetic mechanisms underlie a substantial proportion of variation in each trait. However, no explanatory SNPs were shared between traits, and upper thermal tolerance varied approximately linearly with latitude, whereas hypoxia tolerance exhibited a steep phenotypic break across the contact zone between the subspecies. These results suggest that upper thermal tolerance and hypoxia tolerance are neither phenotypically correlated nor genetically associated, and thus that rates of adaptive change in these traits can be independently fine-tuned by natural selection. This modularity of important traits can underpin the evolvability of organisms to complex future environmental change.

Mitochondria, Temperature, and the Pace of Life.

Life history strategies, physiological traits, and behavior are thought to covary along a "pace of life" axis, with organisms at the fast end of this continuum having higher fecundity, shorter lifespan, and more rapid development, growth, and metabolic rates. Countergradient variation represents a special case of pace of life variation, in which high-latitude organisms occupy the fast end of the continuum relative to low-latitude conspecifics when compared at a common temperature. Here, we use Atlantic killifish (Fundulus heteroclitus) to explore the role of mitochondrial properties as a mechanism underlying countergradient variation, and thus variation in the pace of life. This species is found along the Atlantic coast of North America, through a steep latitudinal thermal gradient. The northern subspecies has faster development, more rapid growth, higher routine metabolic rate, and higher activity than the southern subspecies when compared at a common temperature. The northern subspecies also has greater mitochondrial respiratory capacity in the liver, although these differences are not evident in other tissues. The increased respiratory capacity of liver mitochondria in northern fish is associated with increases in the activity of multiple electron transport complexes, which largely reflects an increase in the amount of inner mitochondrial membrane per mitochondrion in the northern fish. There are also differences in the lipid composition of liver mitochondrial membranes, including differences in cardiolipin species, which could also influence respiratory capacity. These data suggest that variation in mitochondrial properties could, at least in part, underlie variation in the pace of life in Atlantic killifish.

Gene expression plasticity in response to salinity acclimation in threespine stickleback ecotypes from different salinity habitats.

Phenotypic plasticity is thought to facilitate the colonization of novel environments and shape the direction of evolution in colonizing populations. However, the relative prevalence of various predicted patterns of changes in phenotypic plasticity following colonization remains unclear. Here, we use a whole-transcriptome approach to characterize patterns of gene expression plasticity in the gills of a freshwater-adapted and a saltwater-adapted ecotype of threespine stickleback (Gasterosteus aculeatus) exposed to a range of salinities. The response of the gill transcriptome to environmental salinity had a large shared component common to both ecotypes (2159 genes) with significant enrichment of genes involved in transmembrane ion transport and the restructuring of the gill epithelium. This transcriptional response to freshwater acclimation is induced at salinities below two parts per thousand. There was also differentiation in gene expression patterns between ecotypes (2515 genes), particularly in processes important for changes in the gill structure and permeability. Only 508 genes that differed between ecotypes also responded to salinity and no specific processes were enriched among this gene set, and an even smaller number (87 genes) showed evidence of changes in the extent of the response to salinity acclimation between ecotypes. No pattern of relative expression dominated among these genes, suggesting that neither gains nor losses of plasticity dominated the changes in expression patterns between the ecotypes. These data demonstrate that multiple patterns of changes in gene expression plasticity can occur following colonization of novel habitats.

Metabolic and regulatory responses involved in cold acclimation in Atlantic killifish, Fundulus heteroclitus.

Ectotherms often respond to prolonged cold exposure by increasing mitochondrial capacity via elevated mitochondrial volume density [V V(mit,f)]. In fish, higher V V(mit,f) is typically associated with increased expression of nuclear respiratory factor 1 (Nrf1), a transcription factor that induces expression of nuclear-encoded respiratory genes. To examine if nrf1 expression or the expression of other genes that regulate mitochondrial biogenesis contribute to changes in whole-organism metabolic rate during cold acclimation, we examined the time course of changes in the expression of these genes and in metabolic rate in Atlantic killifish, Fundulus heteroclitus. Cold acclimation rapidly decreased metabolic rate, but increased the expression of nrf1 more gradually, with a time course that depended on how rapidly the fish were transitioned to low temperature. Cold-induced nrf1 expression was not associated with increases in biochemical indicators of mitochondrial respiratory capacity, suggesting that cold-induced mitochondrial biogenesis may occur without increases in oxidative capacity in this species. These observations imply that changes in nrf1 expression and metabolic rate due to cold acclimation occur through different physiological mechanisms, and that increases in V V(mit,f) are likely not directly related to changes in metabolic rate with cold acclimation in this species. However, nrf1 expression differed between northern and southern killifish subspecies regardless of acclimation temperature, consistent with observed differences in metabolic rate and V V(mit,f) at 5 °C between these subspecies. Taken together, these results reveal substantial complexity in the regulation of V V(mit,f) and mitochondrial capacity with temperature in fish and the relationship of these parameters to metabolic rate.

Mitochondrial genotype and phenotypic plasticity of gene expression in response to cold acclimation in killifish.

Adjustments of aerobic metabolic processes are critical components of organismal responses to environmental change that require tight coordination between the nuclear and mitochondrial genomes. Intraspecific differences in mitochondrial genotype can affect gene transcription in both genomes. Thus, variation in mitochondrial genotype may be associated with differences in the plasticity of gene expression when organisms are faced with changes in environmental conditions. Cold acclimation is known to result in metabolic responses involving increases in mitochondrial amount and capacity, suggesting that low temperatures may pose a particular challenge when coordinating the functions of the nuclear and mitochondrial genomes. In this study, we utilized RNA-seq to assess transcriptome-wide gene expression in the muscle of Atlantic killifish (Fundulus heteroclitus) from a population that contains segregating variation in mitochondrial genotype. We examined gene expression plasticity in response to 5 °C acclimation and the effects of mitochondrial genotype on this plasticity. Cold acclimation resulted in changes in gene expression consistent with up-regulation of genes involved in many cellular functions, including spliceosomal and proteasomal processes, and with down-regulation of genes involved in extracellular matrix, muscle contraction and oxidative phosphorylation functions. There were few differences in gene expression between killifish with different mitochondrial genotypes: 14 genes demonstrated significant interactions between mitochondrial genotype and acclimation temperature and three genes demonstrated effects of mitochondrial genotype alone. These results indicate that variation in mitochondrial genotype has modest effects on gene expression; the majority of which are revealed as differences in plasticity as a result of environmental change.

Conserved effects of salinity acclimation on thermal tolerance and hsp70 expression in divergent populations of threespine stickleback (Gasterosteus aculeatus).

In natural environments, organisms must cope with complex combinations of abiotic stressors. Here, we use threespine stickleback (Gasterosteus aculeatus) to examine how changes in salinity affect tolerance of high temperatures. Threespine stickleback inhabit a range of environments that vary in both salinity and thermal stability making this species an excellent system for investigating interacting stressors. We examined the effects of environmental salinity on maximum thermal tolerance (CTMax) and 70 kDa heat shock protein (hsp70) gene expression using divergent stickleback ecotypes from marine and freshwater habitats. In both ecotypes, the CTMax of fish acclimated to 20 ppt was significantly higher compared to fish acclimated to 2 ppt. The effect of salinity acclimation on the expression of hsp70-1 and hsp70-2 was similar in both the marine and freshwater stickleback ecotype. There were differences in the expression profiles of hsp70-1 and hsp70-2 during heat shock, with hsp70-2 being induced earlier and to a higher level compared to hsp70-1. These data suggest that the two hsp70 isoforms may have functionally different roles in the heat shock response. Lastly, acute salinity challenge coupled with heat shock revealed that the osmoregulatory demands experienced during the heat shock response have a larger effect on the hsp70 expression profile than does the acclimation salinity.

Warm acclimation improves hypoxia tolerance in Fundulus heteroclitus.

Human activities are increasing both the frequency of hypoxic episodes and the mean temperature of aquatic ecosystems, but few studies have considered the possibility that acclimation to one of these stressors could improve the ability to cope with the other stressor. Here, we used Atlantic killifish, Fundulus heteroclitus, to test this hypothesis. Hypoxia tolerance was measured as time to loss of equilibrium in hypoxia (LOEhyp) at 0.4 kPa oxygen. Time to LOEhyp declined from 73.3 ± 6.9 min at 15 °C to 2.6 ± 3.8 min at 23 °C, and at 30 °C no fish could withstand this level of hypoxia. Prior acclimation to warm temperatures significantly increased time to LOEhyp. Hypoxia tolerance of the southern subspecies of killifish, F. heteroclitus heteroclitus, was greater than that of the northern subspecies, F. heteroclitus macrolepidotus, measured both as critical oxygen tension (Pcrit) and as time to LOEhyp. Warm acclimation offset the negative effects of temperature on time to LOEhyp to a similar extent in the two subspecies. Warm acclimation increased total lamellar surface area of the gill in both subspecies as a result of regression of an interlamellar cell mass (ILCM). However, differences in total lamellar surface area could not explain differences in time to LOEhyp between the subspecies, suggesting that the lower time to LOEhyp of northern fish is related to their higher routine metabolic rate. These data suggest that thermal plasticity in gill morphology can improve the capacity of this species to tolerate hypoxia, and shows how existing plasticity may help organisms to cope with the complex interacting stressors that they will encounter with increasing frequency as our climate changes.

Phenotypic plasticity and divergence in gene expression.

The extent to which phenotypic plasticity, or the ability of a single genotype to produce different phenotypes in different environments, impedes or promotes genetic divergence has been a matter of debate within evolutionary biology for many decades (see, for example, Ghalambor et al. ; Pfennig et al. ). Similarly, the role of evolution in shaping phenotypic plasticity remains poorly understood (Pigliucci ). In this issue of Molecular Ecology, Dayan et al. () provide empirical data relevant to these questions by assessing the extent of plasticity and divergence in the expression levels of 2272 genes in muscle tissue from killifish (genus Fundulus) exposed to different temperatures. F. heteroclitus (Fig. A) and F. grandis are minnows that inhabit estuarine marshes (Fig. B) along the coasts of the Atlantic Ocean and Gulf of Mexico in North America. These habitats undergo large variations in temperature both daily and seasonally, and these fish are known to demonstrate substantial phenotypic plasticity in response to temperature change (e.g. Fangue et al. ). Furthermore, the range of F. heteroclitus spans a large latitudinal gradient of temperatures, such that northern populations experience temperatures that are on average ~10°C colder than do southern populations (Schulte ). By comparing gene expression patterns between populations of these fish from different thermal habitats held in the laboratory at three different temperatures, Dayan et al. () address two important questions regarding the interacting effects of plasticity and evolution: (i) How does phenotypic plasticity affect adaptive divergence? and (ii) How does adaptive divergence affect plasticity?

Thermal acclimation is not necessary to maintain a wide thermal breadth of aerobic scope in the common killifish (Fundulus heteroclitus).

Loss of aerobic scope at high and low temperatures is a physiological mechanism proposed to limit the thermal performance and tolerance of organisms, a theory known as oxygen- and capacity-limited thermal tolerance (OCLTT). Eurythermal organisms maintain aerobic scope over wide ranges of temperatures, but it is unknown whether acclimation is necessary to maintain this breadth. The objective of this study was to examine changes in aerobic scope in Fundulus heteroclitus, a eurythermal fish, after acclimation and acute exposure to temperatures from 5° to 33°C. The range of temperatures over which aerobic scope was nonzero was similar in acclimated and acutely exposed fish, suggesting that acclimation has modest effects on the thermal breadth of aerobic scope. However, in acclimated fish, there was a clear optimum temperature range for aerobic scope between 25° and 30°C, whereas aerobic scope was relatively constant across the entire temperature range with acute temperature exposure. Therefore, the primary effect of acclimation was to increase aerobic scope between 25° and 30°C, which paradoxically resulted in a narrower temperature range of optimal performance in acclimated fish compared to acutely exposed fish. There was only weak evidence for correlations between the thermal optimum of aerobic scope and the thermal optimum of measures of performance (specific growth rate and gonadosomatic index), and indicators of anaerobic metabolism (lactate accumulation and lactate dehydrogenase activity) only increased at high temperatures. Together these data fit many, but not all, of the predictions made by OCLTT.